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, S-wave pairing in neutron matter is studied within an extension of correlated basis function (CBF) 

' theory to include the strong, short range spatial correlations due to realistic nuclear forces and the 

pairing correlations of the Bardeen, Cooper and Schrieffer (BCS) approach. The correlation operator 
contains central as well as tensor components. The correlated BCS scheme of Ref. |1], developed for 
simple scalar correlations, is generalized to this more realistic case. The energy of the correlated pair 
condensed phase of neutron matter is evaluated at the two-body order of the cluster expansion, but 
considering the one-body density and the corresponding energy vertex corrections at the first order 
of the Power Series expansion. Based on these approximations, we have derived a system of Euler 
pH ' equations for the correlation factors and for the BCS amplitudes, resulting in correlated non linear 

gap equations, formally close to the standard BCS ones. These equations have been solved for the 
momentum independent part of several realistic potentials (Reid, Argonne v-n and Argonne ug/) to 
stress the role of the tensor correlations and of the many-body effects. Simple Jastrow correlations 
and/ or the lack of the density corrections enhance the gap with respect to uncorrelated BCS, whereas 
it is reduced according to the strength of the tensor interaction and following the inclusion of many- 
body contributions. 
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k> ■ I. INTRODUCTION 

^. 

5t Superfluidity in neutron matter has been a fascinating topic in many-body physics and astrophysics ever since 

Migdal [2] proposed the possibility of superfluid matter in neutron stars. In the inner crust of the star, ^ 5*0 pairing 
in the low density neutron gas permeating the lattice of neutron rich nuclei may occur and peak at densities much 
lower than the empirical nuclear matter saturation density, po = 0.16 fm~'^. A similar pairing may take place for 
the low concentration proton component in the highly asymmetrical nuclear matter in the star's interior. At higher 
interior densities, neutrons may also pair in the anisotropic ^P2-^F2 partial wave. A realistic evaluation of the 
density regimes where superfluidity takes place and of the strength of the connected energy gaps is needed for a 
quantitative tmderstanding of important features of neutron stars, such as the cooling ratelS,^ and the post-glitch 
relaxation times (3-|6|]. 
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The qualitative aspects of superfluidity were shown to be describable in nuclei f?] and in infinite systems of inter- 
acting fermions [8] by the extension of the theory of superconductivity of Bardeen, Cooper and Schrieffer [9] (BCS). 
In terms of the nucleon-nucleon (NN) interaction, it is the long range attraction of the ^ Sq NN potential that dom- 
inates in the inner crust density regime, allowing for S'-wave pairing. The gap closes with rising density since the 
short range repulsion is more and more effective. Proton superfluidity (or superconductivity) has a similar origin 
in the interior, while higher density ^P2-'^^2 neutron pairing is traced back to non central, tensor and spin-orbit, 
components. In BCS theory Cooper-like pairs allow for superfluidity even in presence of the short range repulsion 
of modern potentials. 

On the other hand, the strong nuclear interaction induces short range correlations in the wave function, which also 
largely screen the repulsion and introduce many-body contributions. These two features have competing effects, 
since the former is expected to increase the gap, whereas the latter may diminish it. Modern many-body theories, 
such as the method of correlated basis functions [10] (CBF), the Bethe-Brueckner-Goldstone expansion IllH (BBG), 
the self-consistent Green's functions theory [12] (SCGF), and lately quantum Monte Carlo [13,] (QMC), can, with 
efficiency and accuracy, deal with short range correlations in normal phase nucleonic matter. It can be reasonably 
expected that these methods also may be able to provide a similarly realistic description of the superfluid phase, 
especially when modern NN potentials are used. 

Within CBF the short range correlations are introduced by acting with a many-body correlation operator on a 
set of model functions, so defining a correlated basis to be used in a perturbative expansion where the highly non 
perturbative short range correlation effects are already embedded in the basis. The zeroth order of the correlated 
perturbative expansion corresponds to a variational approach, since the correlation operator (and the ground-state 
model wave function) can be derived applying the Ritz variational principle. The variational level may already give 
reliable results if the correlation operator is chosen in an appropriate way. Because realistic NN potentials have 
important spin- and isospin-dependent components, both central and non central (e.g the tensor potential, mainly 

originating from one-pion exchange), a good variational choice for the pair correlation f{ij) must include at least 
six components, 

p— 1,6 

where 0(P=i^2,3)(,^-) ^ S{ij) - i3f^{ij)fp{ij)-5,,fj)a^{i)afjij),eindd^P'=P+^^{ij) = 6(f)(ij)®r(i) •T(j). 

The greek indices denote the Cartesian components. This choice of the operatorial dependence of the correlation is 
consistent with the use of the non central and momentum independent vq potentials of the form, 

veiij) = ■ (2) 

p=l,6 

However, fe{ij) is in general a very good variational choice /or all the realistic potentials. The introduction of such 
structures directly in the correlation operators allows the variational approach to describe microscopically the struc- 
ture of nuclear matter 1 14] and finite nuclei 1 15] with a good accuracy. 

In this paper we are only dealing with pure neutron matter (PNM), therefore T{i) ■ T{j) = 1 and the 6-operator 

algebra underlying ve{ij) and feiij) reduces to the first 3 components p — 1, 2, 3, where /p^j^ — /^^^ + and 

Since the operators in ([l]! do not commute, the many-body correlation operator, Fq{1,2, ..N), is given by the 
symmetrized product. 



F6(1,2,..7V) 



n /^(U) 

i<j = l,N 



(3) 



In CBF theory such operators are kept fixed for all the intermediate states. The correlated CBF intermediate states 
are obtained by acting with Fq{1, 2, ..N) on the corresponding uncorrelated Slater determinant. 

An alternative approach, hereafter denoted as CBF-J, consists in starting with a simpler Jastrow correlation (IC 
depending only on the interparticle distance, 

Fj{l,2,..N) = n /^('^y)' (4) 

i<Cj=l.N 

and introducing the spin/isospin dependence via a Jastrow-correlated perturbative expansion fl^. This choice may 
not be very efficient since the whole spin-isospin dependence must be perturbatively included. However, the terms 
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of the CBF-J expansion have a much simpler structure than those of the CBF expansion, based on spin dependent 
correlation operators, and can be computed by Fermi hypernetted chain (FHNC) resummation [Ij. A possible draw- 
back of the CBF-J perturbative expansion is the complexity of going beyond the second order perturbation level 
which may be insufficient in the Jastrow-like CBF theory. 

Variational CBF theory has been applied to the S'-wave nucleonic superfluid in Ref . HtIi using central potentials 
and correlations, without tensor components, 

M^J) = E -^'^\n^)P''^\^j) , (5) 

Uj) = E f^'''Hn,)P'''^H^j) , (6) 

where p(-^^) (ij) are projectors onto the two-body subspace of total spin-isospin ST. The A'^-body correlation oper- 
ator is then given by: 



F4(l,2,..iV) = S 



n 

i<j=l,N 



(7) 



Lowest order cluster expansion was used to derive a correlated gap equation with the V4 version of the Reid soft core 
NN interaction [18, 19]. This correlated theory was developed within the independent Cooper pairs approximation 
and does not consider the dependence of the correlation on the BCS amplitudes. The approach takes essentially 
into account the screening of the core repulsion due to the repulsive part of the correlation, and leads to a larger 
gap than BCS. Chen et al. [16] studied ^-pairing with the Reid vg potential, including the interaction tensor compo- 
nents, using the independent Cooper pairs approximation. They considered a simple Jastrow correlation rather then 
the correlation operator Fq of Eq. [TJ but computed the variational energy at a higher level of the cluster expansion 
through FHNC theory [20] . A reduction of the BCS gap of about 30% was found, attributable, however, to a rather 
poor choice of the Jastrow factor. The authors of Ref. [16] also computed the second order perturbative CBF correc- 
tion to the pairing matrix element on top of the Jastrow estimate. This approach, which should take into account 
medium polarization, led to a dramatic reduction of the gap by ~ 80%, much larger than all the other estimates of 
the polarization effects, and inconsistent with X-ray observations [21]. Inspite of the fact that the matrix elements 
of CBF perturbation theory are easier to compute in a Jastrow correlated basis, its convergence for large non-central 
potentials in such a basis is still to be assessed. 

The independent Cooper pairs approximation was overcome in ref. (l|], hereafter denoted as I, with a Jastrow fully 
correlated BCS theory. In this work we begin to extend the work of I to the case of correlations having spin-isospin 
dependent, with both central and tensor components (/e model). 

The use of a /e correlation does not allow for a complete sum of the FHNC diagrams, very much the same as for 
the case of normal phase. Similarly to that case the massive resummations of diagrams can be performed using the 
single operator chain (SOC) approximation of Ref. [19]. In this paper we limit our attention to study pure neutron 
matter at the two-body level plus vertex corrections of the cluster expansion of {H — jiN), where /i is the chemical 
potential determined by fixing the correct mean value of the particle number operator (or the density, for infinite 

systems) (TV) = X;m(4iam>- 

The one-body density p — {N)/il, and consequently the vertex corrections in (H), will be here computed at the 
first order of the Power Series expansion 1I20I1 . This approximation guarantees in the normal phase the correct density 
normalization, order by order, and introduces a first flavor of many-body effects. The expectation value (H) will 
be computed at the second order of the cluster expansion, which provides a sufficiently good description of the 
short-range correlations. 

Minimization of {H — pN)2 with respect to the correlation functions /g and to the BCS amplitudes leads to a 
coupled set of Euler and gap equations, which we denote as correlated BCS equation. The solution of such equation is 
a preliminary, very important step towards a full calculation, which will include higher order effects in the evalution 
of both (TV) and {H) and second order perturbative corrections following orthogonal CBF theory of ref. [22]. A 
second approach consists in using the Auxiliarly Field Diffusion Monte Carlo (AFDMC) method to calculate the gap 
energy of a finite number of neutrons in the superfluid phase. Such a method has been used to simulate up to 114 
neutrons in a periodical box to evaluate the equation of state at zero temperature in neutron matter in the normal 
phase ^3]- The extension to superfluid phases can be done using the method developed in the recent work jjU] 
in the study of low density Fermi gas in the regime of large scattering length interaction. AFDMC simulations 
of this type crucially depend upon the choice of a guiding function to fix the nodes and the phases of the wave 
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function. Therefore, the BCS amplitudes resulting from solving the correlated BCS equation are a fundamental input 
to the AFDMC simulations. Preliminary results of that simulations performed with 14 neutrons have already been 
published [25]. Besides the derivation of the correlated BCS equation and its solution for several potentials of the vq 
type (like the truncated versions of the Reid [18], Argonne wi4 and Argonne vs [26] potentials) we have evaluated 
the gap energy with and without vertex corrections. The latter to compare with BBG [U] and SCGF lilZI , the former 
to estimate the effect of the three-body terms of which the vertex corrections are the main part. 

The plan of the paper is as follows: in Section 2 the correlated BCS theory for a fe correlation is presented; Section 
3 contains the Euler and correlated gap equations; numerical results and details on the solution of the equations are 
given in Section 4; Section 5 will briefly discuss our results and give conclusions and perspectives 



II. CORRELATED BCS THEORY 



A correlated wave function for the neutron matter superfluid phase is constructed as 

Ivp,) = F |BCS> , (8) 

where the model BCS-state vector is 

|BCS) - Y[{uk + vkalal^^)\0) . (9) 

k 

Uk and Vk are the real variational BCS amplitudes, satisfying the relation u| + i;| = 1, |0) is the vacuum state and aj„ 
is the fermion creation operator in the single-particle state \m ~ k.a) whose wave fimction is 

{x = r,s\al^\Q)^(t)rn=k,<y{x = r,s) = r/„(s)exp(ifc • r) . (10) 

51 is the normalization volume and r/cr=t,i is the spin wave function with spin projection a. The second-quantized 
correlation operator F is written in terms of the N-particle correlation operators, Fn, as 

F = Fm\^T){^T\^ (11) 

where rnjv specifies a set of N single-particle states with N — 0,2,4,... In coordinate representation and for a 
/e-type correlation we have: 

(xi,X2, ..XTvlFArl*^™) = Fq[1,2, ..N) {(j)„i^{xi)(j)m^{x2).4inN{xN)) A ■ (12) 

The suffix A stands for an antisymmetrized product of single-particle wave functions and Fq{\^2, ..N) is the /g 
N-particle correlation operator ((3). 

In I the cluster expansions of the two-body distribution function, g{ri2), 

9{ri2) = ^ (*.K(n)^.,(nX(r2)vI/.,(r2)|vI/,), (13) 

(Tl,(T2 

and of the one-body density matrix, n(rii/), 

"(^"') = 77 E (*.l*i,(n)*.„(riO|*.), (14) 



CTi .a. 



in the Jastrow correlated case were studied. In the above equations, N are normalization constants, [f) and (r) 
are the destruction and creation field operators. 

In I it was proved that g{ri2) and n(rii') are given by the sum of all the linked cluster diagrams, constructed by 
the dynamical correlation lines (hj = fj — 1 for the Jastrow correlation) and the BCS statistical correlations, 

exp{ik ■ r) {k) , (15) 

cxp(zfc • r) u(k)v{k) , (16) 



lv{r) 


V 


r d^k 


Po J 




luir) 


v 


r (fk 


Po J 
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where v is the spin-isospin degeneracy (v ^2 for PNM) and po is the average density of the uncorrelated BCS model, 
given by 

Po = -J^Ak). (17) 

The FHNC equations, derived in I, sum at all orders the cluster diagrams contributing to g{ri2) and n{rui) in the 
Jastrow case. Here we are dealing with a spin-dependent correlation operator of the type Fq of Eq. Q, reduced to 
the PNM case. 

In addition to the complexity introduced by the spin-dependence, the noncommutativity of f{ij) with f{ik) im- 
plies that any given cluster diagram generates as many clusters as the number of possible ordering of the operators 
presented in the diagram. This is a formidable task, which is not been yet solved. Reasonable approximations have 
been devised [19, 27] to sum up the leading cluster terms. Instead of following such schemes we calculate exactly 
the lowest order correlated cluster terms of g(ri2) and n(ri2 = 0). This is justified be the fact that we consider 
short-range correlations and a low density system. Moreover, we are mainly interested to derive the correlated BCS 
equations. 

It is well known that normalization properties are better approximated by the succesive terms of the power series 
expansion [28] namely the expansion in the number of correlation lines. The energy expectation value is instead 
better evaluated using the expansion in the number of particles, or, equivalently, in the density. Such inconsistency 
can be partionaly resolved by performing a full FHNC summation of both quantities in the case if the elementary 
diagrams give the negligible contributions. Here we will calculate n{ri2 = 0) up to the first orger of PS expansion 
and g{r 12) at the two-body cluster level plus the vertex corrections, evaluated at the first order of the PS expansion, 
to be consistent with n(ri2 = 0). 



A. One-body density and vertex corrections 

For a BCS-type trial fxmction the density is given by: 

P-^- n ■ ^^^^ 

Fluctuations with respect to this average vanish in the thermodynamic limit. We stress that the actual density, p, 
differs from po because the correlation operator affects {N) (see I). Therefore, po has to be considered as a variational 
parameter, and p has to be computed self-consistently. 

The calculation of p follows the FHNC scheme of I. We limit our attention to the FHNC diagrams with zero and 
one correlation lines, e.g. those belonging to the first order of the Power Series cluster expansion. Fig. ^ shows the 
first order diagrams. 

The external point, denoted as 1, is represented by an open dot, whereas the internal points are given as black dots. 
The oriented lines represent exchange lu or U, functions, composed as in I, whereis the dashed ones are dynamical 
correlations — 1. Diagram DA has an overall 1/2 symmetry factor, canceling the 2 factor coming from the two 
exchange loops with oposite orientations. 

In the standard FHNC theory for the normal phase, diagrams Dl-DA add up to give zero contribution: Dl is 
canceled by D'5 and D2 by DA. We are left with the uncorrelated zeroth order diagrams, given rise to the Fermi gas 
momentum distribution, Uk = 6(fcF — k), fcp = {QTr"^ p/vY^'^ being the Fermi momentum. The total density correctly 
coincides, at any order of the power series, with that of the uncorrelated Fermi sea, p = pa. 

In correlated BCS theory this cancellation no longer holds, and corrections to the uncorrelated po are found, namely 
P + Po- 

Following the notation of I, diagram Dl is the first order of the vertex correction of the Ud type, and the other three 
diagrams are included into the vertex correction of the C/g type. Keeping only the linear terms, the density p is given 
by 

p^ pi ^ CdPo , (19) 
where Cd is given, in terms of the cluster terms Ud and Ue, by the following equations: 



Ce = l + Ud, 

Cd=Ce + Ue^l + Ud + Ue 



(20) 
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Dl 






FIG. 1: The lowest order diagrams that contribute to the vertex corrections. 



Accordingly, Ce is the correction to vertices which are arrival points of some exchange lines, whereas the complete 
is the correction for vertices only connected to dynamical lines. 

It is straightforward to extend the algebraic methods given in 1 to the case of the spin-dependent correlations. 

We will now discuss the terms associated with the diagrammatic structures (Dl-DA), contributing to Ud and Ug- 
For Dl we get: 



(21) 



The K'^^^ matrix is given by 



100\ /010\ /OOl 

K^i'^ = I 3 , K'^'^ =1-20, K'^^ =001 
006/ \002/ 111 -2 



(22) 



Two terms are associated to £)3: the first term has both ^^-type exchange lines; in the second term one line is of the 
Z„-t5^e. The total contribution is: 



In the normal phase, u{k)v{k) = and ■v'^{k) = Q{kf — k), and Ud + J/eS = 0. 
Similarly to DS, £^2 has two analogous terms: 



where 



U:2 [d^r ll{r) K'^'A'^Pk (f^'\r)f^^Hr) - SaSji) 



(23) 

(24) 

(25) 
(26) 



with A'' = (1,3,6), Pk = (1,1,0) and = (1,-1,0) (note that the factor 1/2 of the spin-exchange operator is 
included in the factor in front of the integrals). 

Diagram D4 has three different exchange patterns giving the contributions: 
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(i) U^4, having all ly-type exchanges, 
with, 

lvv{r) = — [ 7^-4t exp(ifc • r)v^{k) . 
Again, if v^{k) = e{kF - k), then U^^ + U^a = 0- 



(27) 



(28) 



(ii) U^l\ having two lu exchanges joining at the external point 1, while the third exchange line is of the Zi,-type, 



J JUU 



TJV TJVV 



(iii) t/™, having a Z„ exchange joining with a Z„ one at point 1, the third line being of the Z„-type, 

jjuv ^ Pa f ^3 



with. 



The total Ue^ term is the sum 



In conclusion, Ue is given by: 



luv{r) = — I T^-^T exp(ifc • r)u{k)v^{k) 



PoJ (27r)3 



(29) 



(30) 



(31) 



(32) 



(33) 



B. Potential energy 

We perform the calculation of the expectation value of a vq potential at the two-body order for the cluster expan- 
sion, but including also the vertex corrections at the interaction points, 1 and 2. The reason for going beyond the 
simple two-body approximation in the superfluid phase lies in the correlation driven modification of the expectation 
value of the number operator with respect to BCS, as discussed in the previous subsection. 

The vertex corrections lead to a fully factorized term, similiar to that in the Jastrow correlated BCS case of I, plus 
commutator correction terms. 



(Vh _P I 



V \CdJ 



+ ACd + ACe„ + ACe„ , 



where the AC terms are the commutator corrections. 

The Vd and Vev functions coincide with the direct and exchange terms of the normal phase of PW, 



Vd{r) = ^f^'\r)v'^^\r),f^\r)K'^''A'', 



(34) 



(35) 



After performing the spin algebra corresponding to the last term of the first line of (|34] l, we obtain 
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Veu{r) = J2 f^'^{M'Hr)f^^\r)K'^^' K^'^'A'Bi. (36) 

All the vertex structures D1-D4 contribute to Cd, as discussed in the previous subsection. However, only Dl and 
D2 originate commutator contributions. 

The commutator terms AC are calculated following the algebraic methods of ref. [19]. After some lengthy calcu- 
lations we obtain: 



where 
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ACeu = -Ud^J d\i2ll{r 12) 



X K'^'K^'^'^'Biii - 25,a - S,i - Sn)f^'Hri2)v^'Hn2)f^''Hn2) ■ (39) 



PO 



3 



/ dVi3 Y A'f^'Hn,)f<^'Hri3), (40) 

i=2,3 

C^e2 / d'nslUn,) (/W(^13)/'''(n3)-/(''(n3)/<'Hri3)+4/(3)(ri3)/(3)(ri3)) , (41) 



u:2 = / d'n,llin,)f('\n,)f(^\n,) . (42) 



These expressions sum the commutator corrections which are linear in the vertex corrections Dl and D2. The much 
smaller higher order terms have been disregarded. 



C. Kinetic energy 



We will adopt the Jackson-Feenberg (JF) identity to evaluate the kinetic energy [IJ]. The advantage of using this 
form lies in the fact that the JF kinetic energy operator is mainly constructed by the sum of one- and two-body 
operators, the three-body operators being almost negligable. Other forms, like the Pandharipande-Bethe or the 
Clark-Westhaus ones, have large three-body pieces, and need to go beyond our two-body plus vertex corrections 
approximation. 

The kinetic energy expectation value per particle is given by the sum of a one- and a two-body term: 



T1+T2 



where Ti gives the uncorrelated BCS kinetic energy per particle, 

ly f d^k 



2m p J (27r 



I 2 2 

k V 



(k). 



The JF T2 energy is given by: 



T2 



4m 



d\i2 



rpJF 



1 /Ce 
I' \Cd 



Ti'{ri2)li{ri2)-T;!:{r,2)li{r,2) 



1 



{TiJ{r,2)VHl{r,2) - TiJ{r,2)VHl{r,2))] 



ATd + ATev + ATe: 



(43) 



(44) 



(45) 



Td and T^v, corresponding to the direct and exchange terms of the normal case, are: 



rpJF 



(ri2) = (/«(ri2)VV(*H^i2)" (V/«(ri2))') , 

i 

i.k.l 

TiJ{ri2) = J2k^'''A'Pi (/(^'(ri2)/W(ri2) - (5,i4i) ■ (46) 



i,k,l 

Similarly, the w-terms are: 

7;:'/ (^12) ^Y.^'"'^'^' (/^'^(^i2)V'/('H^i2) - V/«(ri2) • V/C^-^lna) 



i,k,l 



The commutators terms are calculated as for the potential: 



AT, ^^{Ud- K-, + C/e"2)P0 / dVi2 (/'^H^i2)VV'^'(n2) - (V/W(n2))') , (48) 

Are,„ = - — C7d^ / dVi2 V if ''^'A'PKS - - Ski - 6n) (49) 

i,k,l 

(/^^>(ri2)VV<'^(ri2) - V/«(ri2) • V/f ^(ri2)) ?^(ri2) - ^/«(ri2)/('=Hri2)V2z2(^^^)' 



2m u J ^-^ 

i.k.l 



Ski - Sn) (50) 



X ((/«(ri2)VV('=Hri2)-V/«(ri2)-V/«(ri2)) ?,!(ri2) - i/P(ri2)/f (ri2)V2/2(,,2) 



We disregard the small three-body contributions to the JF kinetic energy. As for the potential energy, the commu- 
tator terms include only cluster diagrams which are linear in the vertex corrections. 
The energy expectation value, at the two-body order of the cluster expansion, is: 

E2 = ^Ti+T2 + V2. (51) 

{N)i 

III. EULER AND CORRELATED GAP EQUATIONS 

The Euler and the correlated gap equations form a set of coupled equations, whose solution determines the corre- 
lation fimctions and the correlated BCS amplitudes. They result from the variational requirement: 

S,j^., {H~\N) =^ 0. (52) 
In deriving the equations we will use the two-body approximation previously discussed, 

(H) - {H)2 , 

{N) ~ (iV)i = (53) 

where E2 and pi are given in Eqs. ||5T] | and l fT9] |, respectively. 

We will make further approximations, which we believe are accurate enough, but that can be eventually released. 
They consist of: 

(i) neglecting the commutator terms in the derivation of the Euler equation, while keeping them in the calculation 
of E2 and pi, 
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(ii) decoupling the BCS amplitude, v(q), from the correlation functions, /(*^ (r). As a consequence, we neglect the 
implicit dependence on /*^*' (r) in the functional variation with respect to v{q), and viceversa in the derivation 
of the Euler equations for f^^^{r). 

In this way we arrive at an Euler equation of precisely the same algebraic structure as that of the bare BCS scheme, 
with the Hamiltonian containing paired terms only [29]. However, with respect to the ordinary BCS treatment, there 
is the crucial distinction that the pairing force and the single-particle energies are now renormalized by the dynamical 
correlations. Correlations also affect the mean density through the vertex corrections of the BCS/FHNC theory. The 
explicit formulae are given below. 



A. Euler equations for the correlation functions 

Following the PW notation for standard nuclear matter, where the Schrodinger-like equations are written in the 
T, S channels (here we consider the isospin T = \ channel only), the following changes are made with respect to the 
normal phase equations: 

= 0) : in the singlet channel, eq. (3.12) of PW, 



, 2 



<i>s=o,T=i-^/l+(^) {ll{r) + 2ll{r)). (54) 



2 

Ce 



= 1) : in the triplet channel, eq. (3.14) of PW, 

*s=i,T=i ^ (^^ j ll{r). (55) 

The modifications to the spin-orbit equations are not included since we are dealing with a vq model. 



B. Correlated gap equation 

The correlated gap equations is derived from: 

5v(k){pE2 ' pX) = 0. (56) 
The functional variation of the density is given by: 

6v(k)P = {^ + Ud + Ue)5y(k)Po + Po5v(k){Ud + Ue) , (57) 

where 

Sv{k)Po = PoTT^ — i2v{k)6vk) . (58) 
After performing the tedious variations of the vertex terms, Ud and Ue, we arrive at the expression: 

Sy(k)P ^ Eo{k) {6y(^k)Po) , (59) 

where 

Eo{k) = = 1 + C/rf(3 - Av^k)) + [/e3 

Ov{k)P0 

2-v^(k)- 4:V^(k) f (fq , , „, „. 
+ 2uikHk) J (2^^"('^)"(l^ - '^1)^(1^ - 

l-2v'{k) f d^q 3, 
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and Fu{q) is 

Fuiq) = J dVexp(zq • r) ^ K^'^'A^Bi {h{r)fk{r) - 5a5ki) • (61) 
The variation of p plays a role in: 

(i) the term p\ of Eq. ll57t , giving rise to XEo{q) {Sv(k)Po) 

(ii) the direct term of pE2 pEd, where 

pEd^^J d\,2 (Vdin^) - ^Tdin,)^ . (62) 

In this case we get: 2EdEo {Sv(k)Po) 

(iii) the exchange term of pE2 pE^, where 

pE, ^ ~ l-{c,Pof J d\,2{{VeAri2)ll{ri2)~V,^ir,2)ll{n2)) 
-^\Tev{ri2)ll{ri2) - Teu{ri2)ll{ri2) 

-\ {T2.{ri2)VHl{ri2) - T2u{n2)VHl{r,2))\ } ■ (63) 
This variation applies to the Z„(r) and Z„(r) functions appearing in E^, with the result 

5.ik){pE.) - [m - ^(fc) 2^fc)t(fc) } (^^'(^)^o) , (64) 

where 

+ 2^^(crfi?e), (65) 

Ce 

A(A:) =-clj ^u{\k-q\)v{\k-q\)(^l{q)~^{Tl{q) + ^Tl{q)j^ , (66) 

where Vg^(<z) and Vg^(q) are the Fourier transforms of Vev{r) and Veu{r) of Eqs. l|35] l and l|36|l , respectively. 
Similarly, T^^{q) and T^M) are the Fourier transforms of r//(r) and T2'^/(r) of Eq. (|l6j; rX(g) and r2^(g) are 
the Fourier transforms of T;,{f (r) and Tj^f (r) of Eq. (|47ll . 

Notice that S(fc), of Eq. ( |65t , includes the constant term provided by the fimctional variation of the vertex correc- 
tion. 

After collecting all the terms, we may write a correlated gap equation, or Euler equation for the correlated BCS 
amplitudes, in the form: 

XEoik) v{k) + E(fc)w(fc) - A(fc) ^ ^ J= = , (67) 



2to "wy w '2^1-'y2(fc) 



which resembles that obtained in standard BCS theory 11291] . The solution for w^(fc) of the correlated gap equation can 
be written as: 



with 



<k) = ^ + S(A:) - A^o(fc) , (69) 
2m 

E{k) = VA2(fc) + e2(fc) , (70) 

A = A - 2^;^ , (71) 
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where A(fc) has to be interpreted as the correlated gap function. Its value at fc = fcp, Ap, is the energy gap, namely 
the energy required to break a pair at the Fermi surface. In the present case, the functions Eo{k), S(fc) and A(fc) all 
depend upon w^(fc), making the correlated gap equation |(68j highly non linear. 



C. Correlated versus standard BCS equation 

The correlated BCS equation l(68l l has the same algebraic structure as the uncorrelated one (see ref. |l29l], Eqs. (5.29) 
and (5.30)). However, the standard BCS equations do not contain the E(fc) term, whereas in our approach S(fc) ^ 0, 
even if the correlation operator is set equal to 1. In fact, in this case the quantities Eo{k), S(fc) and A(fc) become: 



Eo{k) -> 1 

^(^) ^ - / 7^v''i\k-q\)[vc{q) + iv„{q)] , 



(27r)3^ 

A(fc) - - 1 - QD^dfc - g|) W[q) ■ivM)\ ■ (72) 

In correlated BCS, I](fc) dresses the single particle energies h?k'^/2m, and Eo(k) renormalizes the mass. 
Similarly, A(fc) assumes the role of the gap function. From Eq. (5.32) of ref. [29] 

where 

Vfcq = J d\ exp(z(fc - q) ■ r)V{r) . (74) 
From Eq. l|68|l and the normalization relation, u'^{k) + v'^{k) = 1, it follows that 



Therefore 



A(i:).-5^H,«(5M,), (76) 



coincides with Eq. \ 

The comparison with the uncorrelated BCS theory allows identifying E{k) with the excitation energy of the broken 
pair (BP) with respect to the ground-state, as defined in Ref. [29], 

E^v - Egs = E{k) . (77) 



IV. RESULTS 



We have solved the BCS and correlated BCS equations for neutron matter with a variety of potentials, namely 
the Reid (R), Argonne wi4 (A14) and Argonne wg' (A8') ones. In solving the gap equations we have generalized the 
method described by Khodel et al. in Ref. ||30|[. According to this method the original gap equation is identically 
replaced by a set of coupled equations: a non-singular quasilinear integral equation for the dimensionless profile 
function, xl^)/ defined by A(fc) — x(^) Ap and a non-linear algebraic one for the gap, Ap — A{k-p), at the Fermi 
surface. After integrating Eq. ((66) over the angle , we obtain 



A(fc) = -cl / 



q^dq V{k,q)A{q) 
2vr2 2VA2(9)+e2(q) ' 
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with. 



Vik,q) ^At: I r'drjo{kr){Veu{r) 







Teu{r) ^ T2u{r) 



h^ka 

—^Att r^drji{kr)T2u{r).n{qr). (79) 
2m 







It is assumed that the interaction V{k, q) is different from zero at the Fermi surface, ^(fcF, k-^) ^ 0. To solve the gap 
equation we decompose the potential, V{k, q), into a separable part and a remainder, W{k, q), that vanishes when 
either argument is at the Fermi surface: 

V{k,q) - Vvmm+W{k,q) , (80) 



where M^(fcF, k) = W{k, k^) = and = V{k, k^)/VY- Then, the gap equation (|78l l is readily seen to be equivalent 
to an integral equation for the shape function, x(fc)/ 



m + ci r t^^M^= . m , (81) 



together with the algebraic equation, 

for the gap amplitude Ap (assumed nonzero). Since T4^(fc, fcp) is zero by construction, the integral equation l(8l]l has a 
nonsingular kernel, the log-singularity of the BCS equation having been isolated in the amplitude equation (|82] |. An 
iterative solution of this set of equations converges very rapidly. 

The correlated gap equations are solved using the BCS solution at a given fcp as an input. We find that the final 

density, p, is always very close to the initial one, po • The maximum difference between fc™^"* and /cp is well below 
one percent. In Table (Pi we show the input and output values of fcp, of the density, p, of the chemical potentials, Af, 
of the effective mass, m* /m, defined by the relation: 

-¥[k^l,__, >e(.)~ + E(.), (83) 



m 



and of the gap A p obtained with the A8' model for the uncorrelated BCS, and for the Jastrow (J) and /g correlated 
(CO) cases. 

It is evident that the introduction of the correlations very slightly affects the total density On the contrary the 
chemical potential is reduced by the Jastrow correlations by ^ 20 to ^ 30%. Spin dependent correlations provide a 
further, even if small, decrease of Xp — h^k'^/2m. The effective mass, computed via the self-energy consider- 
ably decreases after the introduction of the correlations. The normal phase effective mass, computed microscopically 
in CBF, at kp = 0.8 fm"^ is - 0.8. 

In Figure l|2) we show the Jastrow, spin and tensor correlations at fcp — 0.6 fm^^ for the A8' potential. The 
dash-dotted lines are the normal phase correlations, whereas the solid lines give the correlations after solving the 
correlated gap equations. For the S'-pairing case Jastrow and tensor correlations do not change from the normal to 
the BCS phases. Instead, the spin correlation shows some sensitivity to the environmental phase. It is reasonable to 
expect that for the ^P2-^F2 pairing the tensor correlation also will depend on the phase. 

The V? , 2uv and amplitudes, both for the pure and correlated BCS cases, are shown in Figure (|3) at three Fermi 
momenta. At the lowest value, fcp = 0.1 fm~^, the uncorrelated and correlated amplitudes substantially differ among 
each other, the correlated ones showing a larger deviation from the step function, consistent with the larger gap value 
(Ag = 0.07 MeV and A^o = o.l4 MeV). At fcp = 0.6 fm"^ the amplitudes are very close in both approaches, yielding 
similar gaps (A^ = 2.27 MeV and A^° = 2.25 MeV). At the largest value, fcp = 0.8 fm"\ the correlated amplitudes 
are practically step functions. In fact this is almost the highest density for which we find solution to the correlated 
BCS equations. 

The gap function, A(fc), at fcp = 0.6 fm^ is given in Figure iSll. In addition to the pure and correlated BCS 
functions, we show the one obtained by a simple Jastrow-correlated wave-function. At low fc-values, k /fcp < 2.5, the 
effects of the Jastrow and spin-dependent correlations compensate, providing a gap function close to the BCS result. 
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At larger momenta they add and the correlated gap ftmction departs from the uncorrelated one, up to k/k-^ ^ 15, 
where all functions have essentially vanished. 

The self-energy, S(fc), is depicted in Figure ISll at fcp 0.6 fm \ The main difference between the BCS and 
correlated BCS cases lies in the sharp rising of the correlated S(fc) at fc fcp, which produces the much lower effective 
mass given m Table ^, {m*/mf = 0.96 and (m*/m)'^° = 0.62. The mass renormalization, caused by short-range 
correlations, enhanced the dispersive effect of the mean field, which leads to quenching of the energy gap, which is 
enhanced by the screening effect of the neutron pairing potential. 

Figure ^ displays the energy gap at the Fermi surface, Ap^-*, as a function of the Fermi momentum for the A8' 
potential in the uncorrelated BCS case. The curves correspond to the full and to the decoupling approximation solu- 
tions of Ref . 1 16, 31], with and without the self-energy insertions of eq. |[72l l. The two gaps are very close for S(fc) = 0, 
whereas, after the introduction of the self-energy, the decoupling approximation appears to slightly overestimate the 
full solution. 

Figure ^ gives the gaps for different types of correlations (Jastrow and /g) and at various levels of the cluster 
expansion for the same potential. The A*' gaps are the standard BCS results, those with the superscript 'J' are 
obtained within the Jastrow correlated theory and the 'CO' superscript denotes the corresponding correlations. The 
2b and 36 subscripts in the correlated gaps refer to the pure two-body cluster case and to the one in which the density 
and the vertex corrections are computed at the first order of the power series expansion of Fig. The inclusion 
of the Jastrow and /g correlations in the 2b case enhance the gap, because the short-range repulsion of the potential 
is renormalized by the short-range correlations. The 36 cases include medium modification effects via higher order 
cluster terms. Their effect is quite sizeable and reverse the behavior, both reducing the density region where we 
find a BCS solution and decreasing the maximum gap with respect to the standard case for the spin-dependent 
correlations. In fact, A°(max) - 2.6 MeV at fcp ~ 0.9 fm"\ while Afjuiax) - 2.2 MeV at fcp - 0.6 fm"\ These 
results indicate that higher order many-body cluster terms may be relevant to estimate the gap. 

Finally, in Figure (|8]l we show the gaps for different potentials in the BCS and /g-correlated theories. We have 
used, besides the Argonne vs' model, also the Reid and Argonne vi4 (A14) potentials. These potentials differ mostly 
for the strength of the one-pion exchange induced components. In fact, A14 has much stronger spin and tensor 
potentials than Reid and A8'. This difference shows up in the gaps, in both approaches. The BCS gap is larger in 
A14 than for the other potentials, and more drastically reduced in the correlated case, where AAi4(max) ~ 1.7 MeV 
at fcp ^ 0.5 fm"\ 

Our results for the A14 potential are qualitatively similar to those of Ref. (32*] (see also Ref. fsS], where the more 
recent calculations was done), where the medium polarization was included via Landau theory. The authors found 
an analogous decrease of the BCS gap, with A^™''^"(max) ~ 1.5 — 2 MeV at /cf ^ 0.8 fm^^, but with a wider density 
region allowing for a superfluid solution. 

V. CONCLUSIONS AND PERSPECTIVES 

The problem of an accurate determination of the BCS gap in a strongly interacting matter of nucleons is a long- 
standing one. Medium modification effects are expected to be important, but of difficult quantitative evaluation. 
We have used FHNC/BCS theory to take care of the short range correlations induced by the interaction in neutron 
matter at zero temperature. We have adopted the realistic Argonne I's' two-nucleon potential and a correlation 
factor having central, spin and tensor dependent components. The density has been computed at the first order 
of the power series expansion, since this expansion provides at each order the correct density normalization in the 
normal phase. Consistently, the matrix elements of the hamiltonian in the correlated BCS state are evaluated at the 
two-body cluster level plus vertex corrections at the interacting pair. This treatment, in conjunction with the use of 
spin and tensor correlations lowers the maximum gap at fcp by ^ 20% with respect to the uncorrelated BCS case. 
Moreover, Ap is shifted to a lower density. It is clear from our results the relevance of state dependent correlations 
for a reliable estimate of Ap in neutron matter, as well as the need for inserting medium modifications via higher 
order terms of the cluster expansion. Simple Jastrow, spin independent correlations always enhance Ap, even if mas- 
sive summations of cluster diagrams are performed. This effect is due to the screening of the short-range repulsive 
interaction provided by the Jastrow correlations. Similar conclusions are drawn when the short-range correlations 
are introduced by medium effects within the Brueckner G-matrix theory. State dependent correlations reverse this 
scenario and, after the inclusion of the vertex corrections, reduce Ay. A qualitatively analogous result is found when 
state dependence is introduced by a CBF based perturbative expansion theory on top of Jastrow correlated states, 
but with spin dependent interactions. 

In conclusion, we have stressed in this paper the importance of state dependent correlations and medium effects 
in superfluid neutron matter. Both of these tend to reduce the ^5*0 pairing gap, confirming previous studies. 
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TABLE I: Fermi momentum, fcp, density, p, chemical potential, Ap, effective mas, ra* jra, and gap value, Ap, in different approxi- 
mations (see text). 



The large effects found either by extending the study to the full FHNC/SOC calculations or after the introduction 
of the vertex corrections, strongly point to the need of a realistic estimate of many body effects. This can be done 
by using the calculated correlated BCS amplitudes as the guiding function of an AFDMC calculation. In the latter 
approach it is crucial to have a realistic guiding function in the path constraint. An extension of AFDMC to deal with 
the superfluid phases of neutron matter has been recently made and preliminary results, obtained for 14 neutrons, 
are given in Ref. (l^. A full description of both the AFDMC/BCS method and the corresponding results obtained for 
large systems will be given in a forthcoming paper [34]. An important issue is the role of the long range correlations. 
This can be most easily done in a Jastrow correlated BCS case. Work in this direction is in progress. 
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FIG. 2: (colored online) The central, spin and tensor correlation functions at the fcp = 0.6 fm^^ for the A8' potential. The dash- 
dotted, blue lines are the normal phase correlations. The dashed and solid, black lines are the BCS correlations without and with 
vertex corrections, correspondingly. 
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FIG. 3: (colored online) v'^{k) (solid, black lines), 2v{k)u{k) (dash-dotted, red lines) and it^(fc) (dashed, blue lines) amplitudes 
obtained from the A8' potential at three densities. The thin lines are the pure BCS results; the thick ones represent the correlated 
BCS amplitudes. 
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FIG. 5: (colored online) The self-energy E(A:), in MeV, for the BCS (dash-dotted, blue line) and correlated BCS (solid, black line) 
cases for the A8' model at fcp = 0.6 fm^^. 
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FIG. 6: ^So pairing gaps obtained from the A8' potential in the uncorrelated BCS case. The curves labeled ^deco ^^'^ ^dec 
calculated in the decoupling approximation, without and with self-energy insertions, respectively. The remaining curves refer to 
the full solution (Aj! with E(A:) = and A" with E{k) / 0). 
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FIG. 7: A8' ^So pairing gaps for different correlations and levels of the cluster expansions. See text. 
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FIG. 8: ^Sq pairing gaps for different nucleon-nucleon potentials. 



